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QUARKONIUM PHYSICS AND EFTS 

i ■ 

Effective field theories for the description of quarkonium processes have been newly 
developed and are being developed and provide a unifying description as well as a solid 
and versatile tool giving well definite predictions [QIM 0]. They rely on one hand on 
t-h ■ high order perturbative calculations and on the other hand on lattice simulations, the 

recent progress in both fields having added a lot to the theory reach. Heavy quarkonium 
is a multiscale system and as such is on one hand particularly suitable to be treated in 
an effective field theory framework. On the other hand the existence of many scales 
in quarkonium makes it a unique system to study complex environments. Quarkonium 
probes all the regimes of QCD, from the high energy region, where an expansion in the 
the coupling constant is possible, to the low energy region, where nonperturbative effects 
dominate. It probes also the intermediate region between the two regimes. In particular 
for quarkonium system with a very small radius the interaction turns out to be purely 
perturbative while for system with a large radius with respect to the confinement scale 
the interaction turns out to be nonperturbative. In the complex environment of heavy 
ion collisions quarkonium suppression constitutes a unique probe of deconfinement and 
quark gluon plasma formation |4|]. The different radius of the different quarkonia states 
induces the phenomenon of sequential suppression, allowing to use quarkonium as a 
kind of thermometer for the measurement of the temperature of the formed medium. On 
similar ground quarkonium may constitute a special probe to be used in the study of a 
nuclear medium. The large mass, the clean and known decays mode make quarkonium 
an ideal probe of new physics in some well defined window of parameters of physics 
beyond the Standard Model |0,H. 

The modern approach to quarkonium physics consists in taking advantage of the 
hierarchy of non-relativistic energy scales in the system by constructing a suitable 
hierarchy of effective field theories (EFTs) [Q]]. 

The energy scales are: the heavy-quark mass (hard scale), m, the typical momentum 
transfer (soft scale), p ~ mv, whose inverse sets the typical distance, r, between the 
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heavy quark and the antiquark, and the typical kinetic energy (ultrasoft scale), E ~ mv 2 , 
whose inverse sets the typical time scale of the bound state. The heavy-quark bound- 
state velocity v is a small quantity v 1 (v 2 ~ 0.1 for bb, v 2 ~ 0.3 for cc, v 2 ~ 0.01 
for it), the mass is a large quantity m ^> Aqcd> 0C s (m) <C 1. For energy scales close to 
Aqcd> perturbation theory breaks down and one has to rely on nonperturbative methods. 
Regardless of this, the nonrelativistic hierarchy of scales: m>p~ l/r~mv>£~ mv 2 
also persists below the Aqcd threshold. While the hard scale is always larger than Aqcd> 
different situations may arise for the other two scales depending on the considered 
quarkonium system. The soft scale, proportional to the inverse typical radius r, may 
be a perturbative Aqcd) or a nonperturbative scale (~ Aqcd) depending on the 
physical system. The first case is likely to happen only for the lowest charmonium and 
bottomonium states. We do not have direct information on the radius of the quarkonia 
systems, and thus the attribution of some of the lowest bottomonia and charmonia states 
to the perturbative or the nonperturbative soft regime is at the moment still ambiguous. 
The ultrasoft scale may still be perturbative only in the case of tt threshold states. 
All quarkonium scales get entangled in a typical amplitude involving a quarkonium 
observable. In particular, quarkonium annihilation and production happen at the scale m, 
quarkonium binding happens at the scale mv, which is the typical momentum exchanged 
inside the bound state, while very low-energy gluons and light quarks (also called 
ultrasoft degrees of freedom) live long enough that a bound state has time to form and, 
therefore, are sensitive to the scale mv 2 . Ultrasoft gluons are responsible for phenomena 
like the Lamb shift in QCD. 

A hierarchy of EFTs may be constructed by systematically integrating out modes 
associated to high energy scales not relevant for quarkonium. Such integration is made 
in a matching procedure enforcing the equivalence between QCD and the EFT at a 
given order of the expansion in v. The EFT realizes a factorization at the Lagrangian 
level between the high energy contributions, encoded into the matching coefficients, and 
the low energy contributions, carried by the dynamical degrees of freedom. Poincare 
symmetry remains intact in a nonlinear realization at the level of the NR EFT and 
imposes exact relations among the matching coefficients [8|]. 

At the scale m the suitable EFT is NRQCD H], which follows from QCD by in- 
tegrating out the scale m. As a consequence, the effective Lagrangian is organized as 
an expansion in 1/m and 0C s (m). The field of quarkonium production has seen terrific 
progress in the last few years both in theory and in experiments, for a review see |0, l- 

For what concerns decays, recently, substantial progress has been made in the evalua- 
tion of the factorization formula at order v 7 [SO, in the lattice evaluation of the NRQCD 
matrix elements ffioh and in the data of many hadronic and electromagnetic decays Jj]. 
The data are clearly sensitive to NLO corrections in the Wilson coefficients and pre- 
sumably also to relativistic corrections. Improved theory predictability would entail the 
lattice calculation or data extraction of the NRQCD matrix elements and perturbative 
resummation of large contribution in the NRQCD matching coefficients. The new data 
on hadronic transitions and hadronic decays pose interesting challenging to the theory. 

Lattice NRQCD calculations have undergone a steady development in last few years 
see [0,1. 



At the scales mv and mv 2 the suitable EFT is potential NonRelativistic QCD (pN- 
RQCD) BatZD, which follows from NRQCD by integrating out the scale mv. As a con- 
sequence, the effective Lagrangian is organized as an expansion in 1/m and a s (m), in- 
herited from NRQCD, and an expansion in r 

^pNRQCD = / j3 r ££ C n(«s(M),M) Vn)fe(r;iU>) yk x ^(^V 2 ,...), 
J n k 

where O k are the operators of pNRQCD that live at the low-energy scale mv 2 , /! 
is the NRQCD factorization scale, fl' is the pNRQCD factorization scale and V n k 
are the Wilson coefficients of the EFT that encode the contributions from the scale r 
and are non-analytic in r. Looking at the equations of motion of pNRQCD, we may 
identify V^o with the 1/m" potentials that enter the Schrodinger equation and V n ^o 
with the couplings of the ultrasoft degrees of freedom, which provide corrections to the 
Schrodinger equation. This EFT is close to a Schrodinger-like description of the bound 
state. The bulk of the interaction is carried by potential-like terms V n ,o, but non-potential 
interactions V n ^Q, associated with the propagation of low-energy degrees of freedom 
(QQ colour singlets, QQ colour octets and low energy gluons), are generally present. 
They start to contribute at NLO in the multipole expansion of the gluon fields and are 
typically related to nonperturbative effects ||7D. 

In what follows we will focus on the EFT at the scale mv and mv 2 . Then, there are 
several cases for the physics at hand. In the case in which the EFT has been constructed 
dzl (3 [fl], i.e. for states below threshold, the work is currently going in calculating 
higher order perturbative corrections in v and a s for processes of interest, resumming 
the logarithms in the ratio of the scales that may be sizeable, calculating or extracting 
nonperturbatively low energy correlators and extending the theory with the addictions 
of electromagnetic effects fl2lll and the consideration of QQQ and QQq systems lfl5ll . 
The issue here is precision physics and the study of confinement. Close to threshold the 
EFT has not yet been constructed and the degrees of freedom have still to be identified 
[|3ll [3|l. At finite temperature the EFT is being constructed and the existing results hint 
at a new physical picture with possible application at heavy ion collisions at LHC. 

Below we will review new results focusing on the example of the calculation of the 
interquark potential at zero and at finite temperature in the different dynamical situations. 

QUARKONIUM POTENTIAL AT ZERO TEMPERATURE 

For states away from threshold we have a clear effective field description called pN- 
RQCD, based on perturbative and lattice computations. This is nowadays the standard 
description. 

The soft scale mv may be larger or not than the confinement scale Aqcd depending 
on the radius of the quarkonium system. When mv 2 ~ Aqcd> we speak about weakly- 
coupled pNRQCD because the soft scale is perturbative and the matching from NRQCD 
to pNRQCD may be performed in perturbation theory. When mv ~ Aqcd> we speak 
about strongly-coupled pNRQCD because the soft scale is nonperturbative and the 
matching from NRQCD to pNRQCD may not be performed in perturbation theory. 



The potential is a Wilson coefficient of an EFT. In general undergoes renormalization, 
develops scale dependence and satisfies renormalization group equations which allow to 
resum large logarithms. 



weakly-coupled pNRQCD 

If the quarkonium system is small, the soft scale is perturbative and the potentials can 
be entirely calculated in perturbation theory [HI]. 

Since the degrees of freedom that enter the Schrodinger description are in this case 
both QQ color singlet and QQ color octets, both singlet and octet potentials exist. 
The static singlet QQ potential is pretty well known. The three-loop correction to the 
static potential is now completely known: the fermionic contributions to the three-loop 
coefficient [12211 first became available, and more recently the remaining purely gluonic 
term has been obtained ||23L l24fl . 



The first log related to ultrasoft effects arises at three loops ft 1 311 . Such logarithm 
contribution at N 3 LO and the single logarithm contribution at N 4 LO may be extracted 
respectively from a one-loop and two-loop calculation in the EFT and have been calcu- 
lated in 03 • 

The perturbative series of the static potential suffers from a renormalon ambiguity 
(i.e. large /3o contributions) and from large logarithmic contributions. The singlet static 
energy, given by the sum of a constant, the static potential and the ultrasoft corrections, 
is free from ambiguities of the perturbative series. By resumming the large logs using 
the renormalization group equations and comparing it (at the NNLL) with lattice cal- 
culations of the static energy one sees that the QCD perturbative series converges very 
nicely to and agrees with the lattice result in the short range (up to 0.25 fm) and that no 
nonperturbative linear ("stringy") contribution to the static potential exist [|14l ll2ll. 

In particular, the recently obtained theoretical expression for the complete QCD 
static energy at NNNLL precision has been used to determine roA^ s by comparison with 
available lattice data, where tq is the lattice scale and Ajj s is the QCD scale, obtaining 
roA^ 5 = 0.622+QQjg for the zero-flavor case. This extraction was previously performed 
at the NNLO level (including an estimate at NNNLO) in tf25|] . The same procedure can 
be used to obtain a precise evaluation of the unquenched roAjj S value after short distance 
unquenched lattice data for the QQ exist. 

The static octet potential is known up to two loops rf26ll . Relativistic corrections to the 
static singlet potential have been calculated over the years and are summarized in [Q]]. 

In the case of QQq baryons, the static potential has been determined up to NNLO 
in perturbation theory |fl5|l and recently also on the lattice [|29ll . Terms suppressed by 
powers of 1/m and r in the Lagrangian have been matched (mostly) at leading order and 
used to determine, for instance, the expected hyperfine splitting of the ground state of 
these systems. 

In the case of QQQ baryons, the static potential has been determined up to NNLO 
in perturbation theory [|15|l and also on the lattice [|28|]. The transition region from a 
Coulomb to a linearly raising potential is characterized in this case also by the emergence 
of a three-body potential apparently parameterized by only one length. It has been shown 



that in perturbation theory a smooth genuine three-body potential shows up at two loops. 



strongly-coupled pNRQCD 

If the quarkonium system is large, the soft scale is nonperturbative and the potentials 
cannot be entirely calculated in perturbation theory yj]. Then the potential matching 
coefficients are obtained in the form of expectation values of gauge-invariant Wilson- 
loop operators. In this case, heavy-light meson pairs and heavy hybrids develop a mass 
gap of order Aqcd with respect to the energy of the QQ pair, the second circumstance 
being apparent from lattice simulations. Thus, away from threshold, the quarkonium 
singlet field S is the only low-energy dynamical degree of freedom in the pNRQCD 
Lagrangian [fl6L [fl, [l8|] (neglecting ultrasoft corrections coming from pions and other 
Goldstone bosons). The singlet potential Vs(r) can be expanded in powers of the inverse 
of the quark mass; static, 1/m and 1/m 2 terms were calculated long ago [16]. They in- 
volve NRQCD matching coefficients (containing the contribution from the hard scale) 
and low-energy nonperturbative parts given in terms of static Wilson loops and field- 
strength insertions in the static Wilson loop (containing the contribution from the soft 
scale). Such expressions correct and generalize previous finding in the Wilson loop ap- 
proach II 1711 that were typically missing the high energy parts of the potentials, encoded 
into the NRQCD matching coefficients and containing the dependence on the logarithms 
of m, and some of the low energy contributions. Poincare invariance [8] establishes ex- 
act relations between the potentials of the type of the Gromes relation between spin 
dependent and static potentials [|19|l . 

In this regime of pNRQCD, we recover the quark potential singlet model. However, 
here the potentials are calculated in QCD by nonperturbative matching. Their evaluation 
requires calculations on the lattice or in QCD vacuum models [13211 . 

Recent progress includes new, precise lattice calculations of these potentials obtained 
using the Liischer multi-level algorithm Il27ll . 

The nonperturbative potentials for the QQQ and QQq have been obtained in the 
second reference of IU5I1 and in (13011 . 



QUARKONIUM CLOSE OR ABOVE THRESHOLD 

In the most interesting region, the region close to threshold where many new states, 
conceivably of an exotic nature have been recently discovered, no EFT description has 
yet been constructed nor the appropriate degrees of freedom clearly identified lt3ll HI]. 
An exception is constituted by the X (3872) that displays universal characteristics related 
to its being so close to threshold, reason for which a beautiful EFT description could be 



obtained [20] 



The threshold region remains troublesome also for the lattice, althought several ex- 
cited states calculations have been recently being pionereed. 

Lattice results about the the crosstalk of the static potential with a pair of heavy-light 



mesons in the lattice have recently appeared fl3 3h but further investigations appear to be 
necessary. 



QUARKONIUM POTENTIAL AT FINITE TEMPERATURE 



The study of quarkonium in media has recently undergone crucial developments. Large 
datasets from heavy-ion collisions have recently become available at RHIC displaying 
new features related to the quark gluon plasma formation characteristics like the particu- 
lar structure of jet quenching and the very low viscosity to entropy ratio. In particular the 
quark gluon plasma looks more like a liquid than a plasma and the use of perturbative 
expansion appears to be justified only at temperature bigger than the deconfinement one. 

The suppression of quarkonium production in the hot medium remains one of the 
cleanest and most relevant probe of deconfined matter. 

However, the use of quarkonium yields as a hot-medium diagnostic tool has turned 
out to be quite challenging for several reasons. Quarkonium production has already 
been found to be suppressed in proton-nucleus collisions by cold-nuclear-matter effects, 
which themselves require dedicated experimental and theoretical attention. Recombi- 
nation effects may play an additional role and thus transport properties may become 
relevant to be considered. Finally, the heavy quark-antiquark interaction at finite tem- 
perature T has to be obtained from QCD. 

For observables only sensitive to gluons and light quarks, a very successfull EFT 
called Hard Thermal Loop (HTL) effective theory has been derived in the past [|430 
by integrating out the hardest momenta propotional to T from the dynamics. However, 
considering also heavy quarkonium in the hot QCD medium, one has to consider in 
addition to the thermodynamical scales in T also the scales of the nonrelativistic bound 
state and the situation becomes more complicate. 

In the last few years years, there has been a remarkable progress in constructing EFTs 
for quarkonium at finite temperature and in rigorously defining the quarkonium poten- 
tial. In [33,[35|], the static potential was calculated in the regime T ^> 1/r > nip, where 
mj) is the Debye mass and r the quark-antiquark distance, by performing an analyti- 
cal continuation of the Euclidean Wilson loop to real time. The calculation was done 
in the weak-coupling resummed perturbation theory. The imaginary part of the gluon 
self energy gives an imaginary part to the static potential and hence a thermal width to 
the quark-antiquark bound state. In the same framework, the dilepton production rate 
for charmonium and bottomonium was calculated in (3a In [|38|l . static particles 
in real-time formalism were considered and the potential for distances l/r ~ mo was 
derived for a hot QED plasma. The real part of the static potential was found to agree 
with the singlet free energy and the damping factor with the one found in [f34^ | . In [39], 
a study of bound states in a hot QED plasma was performed in a non-relativistic EFT 
framework. In particular, the hydrogen atom was studied for temperatures ranging from 
T <C ma 2 to T ~ m, where the imaginary part of the potential becomes larger than the 
real part and the hydrogen ceases to exist. The same study has been extended to muonic 
hydrogen in ll40ll . providing a method to estimate the effects of a finite charm quark mass 
on the dissociation temperature of bottomonium. 

An EFT framework in real time and weak coupling for quarkonium at finite tempera- 
ture was developed in (4211 working in real time and in the regime of small coupling g, 
so that gT <^T and v ~ (X s , which is expected to be valid for tightly bound states: T(15), 
J/\jf, ... . 

Quarkonium in a medium is characterized by different energy and momentum scales; 



there are the scales of the non-relativistic bound state that we have discussed at the 
beginning, and there are the thermodynamical scales: the temperature T, the inverse 
of the screening length of the chromoelectric interactions, i.e. the Debye mass mo and 
lower scales, which we will neglect in the following. 

If these scales are hierarchically ordered, then we may expand physical observables 
in the ratio of such scales. If we separate explicitly the contributions from the different 
scales at the Lagrangian level this amounts to substituting QCD with a hierarchy of 
EFTs, which are equivalent to QCD order by order in the expansion parameters. As 
it has been described in the previous sections at zero temperature the EFTs that follow 
from QCD by integrating out the scales m and mv are called respectively Non-relativistic 
QCD (NRQCD) and potential NRQCD (pNRQCD). We assume that the temperature 
is high enough that T 3> gT ~ mo holds but also that it is low enough for T <C m 
and 1 jr ~ mv > mo to be satisfied, because for higher temperature the bound state 
ceases to exist. Under these conditions some possibilities are in order. If T is the next 
relevant scale after m, then integrating out T from NRQCD leads to an EFT that we may 



name NRQCDhtl* because it contains the hard thermal loop (HTL) Lagrangian 114311 . 
Subsequently integrating out the scale mv from NRQCDhtl leads to a thermal version 
of pNRQCD that we may call pNRQCDhtl- If me next relevant scale after m is mv, 
then integrating out mv from NRQCD leads to pNRQCD. If the temperature is larger 
than mv 2 , then the temperature may be integrated out from pNRQCD leading to a new 
version of pNRQCDhtl [44] . Note that, as long as the temperature is smaller than the 
scale being integrated out, the matching leading to the EFT may be performed putting 
the temperature to zero. 

The derived potential V describes the real-time evolution of a quarkonium state 
in a thermal medium. At leading order, the evolution is governed by a Schrodinger 
equation. In an EFT framework, the potential follows naturally from integrating out all 
contributions coming from modes with energy and momentum larger than the binding 
energy. For T < V the potential is simply the Coulomb potential. Thermal corrections 
affect the energy and induce a thermal width to the quarkonium state; these may be 
relevant to describe the in medium modifications of quarkonium at low temperatures. 
For T > V the potential gets thermal contributions, which are both real and imaginary. 

General findings in this picture are: 

• The thermal part of the potential has a real and an imaginary part. The imaginary 
part of the potential smears out the bound state peaks of the quarkonium spectral 
function, leading to their dissolution prior to the onset of Debye screening in 
the real part of the potential (see, e.g. the discussion in |41|]). So quarkonium 
dissociation appears to be a consequence of the appearance of a thermal decay 
width rather than being due to the color screening of the real part of the potential; 
this follows from the observation that the thermal decay width becomes as large as 
the binding energy at a temperature at which color screening may not yet have set 
in. 

• Two mechanisms contribute to the thermal decay width: the imaginary part of the 
gluon self energy induced by the Landau-damping phenomenon (existing also in 
QED) [|34|l and the quark-antiquark color singlet to color octet thermal break up (a 
new effect, specific of QCD) |46|]. Parametrically, the first mechanism dominates 



for temperatures such that the Debye mass mo is larger than the binding energy, 
while the latter dominates for temperatures such that mo is smaller than the binding 
energy. 

The obtained singlet thermal potential, V, is neither the color-singlet quark- 
antiquark free energy nor the internal energy. It has an imaginary part and may 
contain divergences that eventually cancel in physical observables Q460 . 
Temperature effects can be other than screening, typically they may appear as 
power law corrections or a logarithmic dependence Il46l l39ll. 

The dissociation temperature goes parametrically as 7rr me i t i ng ~ mg^ rf39l l4lll. 



The EFT provides a clear definition of the potential and a coherent and systematical 
setup to calculate masses and widths of quarkonium at finite temperature. In l45ll heavy 
quarkonium energy levels and decay widths in a quark-gluon plasma, below the melting 
temperature at a temperature T and screening mass mo satisfying the hierarchy m(X s 3> 
kT ^> mC( s 2 ^> mo, have been calculated at order ma s 5 . This situation is relevant for 
bottomonium 15 states (T(15), rjb) at the LHC. It has been found Il45ll that: at leading 
order the quarkonium masses increase quadratically with T which in turn implies the 
same functional increase in the energy of the dileptons produced in the electromagnetic 
decays; a thermal correction proportial to T 2 appears in the electromagnetic quarkonium 
decay rates; at leading order a decay width linear with the temperature is developed 
which implies a tendency to dissolve by decaying to the continuum of the colour-octet 
states. 

In [|46l . l47ll the Polyakov loop and the correlator of two Polyakov loops at finite 
temperature has ben calculated at next-to-next-to-leading order in the weak coupling 
regime and at quark-antiquark distances shorter than the inverse of the temperature and 
for Debye mass larger than the Coulomb potential. The calculation has been performed 
also the in EFT framework Ii46ll and a relation between the Polyakov loop correlator and 



the singlet and octet quark-antiquark correlator has been established in this setup [42]. 

First attempts to generalize this new picture to the nonperturbative regime have been 
undertaken in |48|]. 
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